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1.0  INTRODUCTION 

Development  of  numerical  wave  prediction 
models  for  purposes  of  wave  forecasting  and 
hindcasting  has  been  a  key  part  of  wave  research  for 
several  decades.  Models  generally  address  particular 
wave  processes  such  as  wave  generation  and 
propagation  (in  deep  and/or  shallow  water),  wave 
refraction/diffraction,  or  wave  breaking.  Each  of  these 
processes  involves  different  physics,  spatial  scales,  and 
numerical  approaches.  New  types  of  amphibious 
systems  and  strategies  require  an  integrated  suite  of 
models  that  provide  predictive  capability  over  a  large 
region  from  deep-water  to  the  beach  and  along  the 
coast.  Several  state-of-the-art  models  have  been 
developed  to  the  point  that  they  can  be  implemented 
operationally  or  are  already  used  operationally. 

The  Surf  Zone  Environmental  Representations 
(SZER)  Program  is  identifying,  linking,  and  operating 
a  coupled  suite  of  wave  and  surf  models  to  provide 
automated  calculations  of  wave  conditions  from  deep¬ 
water  to  and  along  the  beach.  Unlike  traditional  wave 
forecasting  and  hindcasting,  this  effort’s  goal  is  to 
develop  a  methodology  so  that  wave  conditions  can  be 
calculated  realistically  over  large  regions  for 
simulations  of  military  systems  and  amphibious 
operations.  Outputs  ultimately  will  be  used  as  inputs  to 
simulation  and  visualization  software.  The  developed 
methodology  accommodates  replacement  of  particular 
models  with  updated  or  improved  models  as  they 
become  available.  Models  that  are  incorporated  into 
SZER  are  WAM,  STWAVE,  REFDIFl,  and  the  Navy 
Standard  Surf  Model  (NSSM).  Model  selection  based 
physics  and  practical  implementation 
^considerations,  regions  in  which  the  models  are  most 
applicable,  and  model  coupling  techniques  are 
described.  Examples  of  coupled  model  output  for 


Onslow  Bay,  NC,  an  important  area  for  the  military 
training  and  exercise  community,  are  discussed. 

2,0  MODEL  DESCRIPTION 

The  SZER  approach  is  to  integrate  a  suite  of 
physically  consistent  wave  and  surf  models  originating 
in  deep-water  and  progressing  to  shallow  water  into  the 
surf  zone.  Figure  1  illustrates  the  SZER  approach 
implemented  for  Onslow  Bay,  NC  in  which  a  regional 
WAM  model  is  independently...^  coupled  to  the 
STWAVE  and  REFDIFl  shallow-water  wave  models. 
Outputs  from  one  or  both  of  these  models  are  used  as 
inputs  to  NSSM.  In  this  paper,  comparisons  between 
STWAVE  and  REFDIFl,  as  well  as  recommendations 
for  their  implementation  are  made. 


Figure  1  SZER  modeling  approach  in  which  the  deepwater 
WAM  model  is  coupled  independently  to  the  STWAVE  and 
REFDIF  shallow  water  wave  models;  the  outputs  from  one  or 
both  of  these  models  arc  used  as  inputs  into  NSSM. 

2.1  WAM 


The  WAM  wave  model  is  a  spectral  wave 
prediction  model  developed  by  the  WAMDI  Group 
(1988;  also  Komen  et  al.  1994),  an  international 
consortium  of  wave  modelers.  WAM  describes  the  sea 
surface  as  a  discretized  two-dimensional  (2-D) 
spectrum  of  sea  surface  elevation  variance  density. 


The  Fleet  Numerical  Meteorology  and  Oceanography 
Command  and  Naval  Oceanographic  Office 
(NAVOCEANO)  run  operational  global  and  regional 
implementations  of  WAM  Cycle  4  (Wittmann  and 
Farrar  1 997).  WAM  used  in  this  study  is  the  Carolina 
regional  WAM  run  by  NAVOCEANO  with  a  0.2® 
resolution  and  spectra  saved  at  selected  locations. 
This  regional  WAM  is  coupled  to  the  1®  North  Atlantic 
WAM.  This  implementation  of  WAM  is  not  run 
routinely  by  NAVOCEANO.  Rather,  it  is  used  for 
exercise  support. 

WAM  is  discretized  into  25  frequency  bands 
with  center  frequencies  ranging  from  0.04333  Hz  to 
0.32832  Hz,  with  each  frequency  being  1.1  times  that 
of  the  next  lower  band.  Direction  is  discretized  into  24 
bands  of  width  15°.  WAM  computes  the  wind¬ 
generated  energy  density  of  each  spectral  component. 
Energy  is  also  propagated  in  space,  with  refraction  due 
to  depth  variation,  and  dispersion  due  to  the  nature  of 
the  waves.  Because  WAM  spatial  resolution  does  not 
resolve  bathymetric  variations  close  to  a  coast,  WAM’s 
refraction  calculations  apply  to  offshore  regions  rather 
than  to  the  regions  covered  by  STWAVE  and 
REFDIFI. 

In  this  study,  a  1 0-day  period  ranging  from  March 
12-22,  1997  was  selected  in  which  WAM  wave  spectra 
were  saved  and  subsequently  used  as  inputs  to  the 
STWAVE  and  REFDIFI  shallow- water  wave  models. 
Figure  2  shows  the-W.^M  domain  and  locations  of  the 
spectra  to  be  used  as  possible  inputs  to  STWAVE  and 
REFDIFI. 
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2.1.1  WAM  Wind  Forcing 

The  Navy  Operational  Regional  Atmospheric 
Prediction  System  (NORAPS)  for  the  Continental 
United  States  (CONUS)  10-m  winds  were  used  as 
inputs  into  the  WAM  (and  STWAVE)  hindcasts 
performed  for  the  area  including  Onslow  Bay,  NC 
COlWS  NORAPS  has  a  horizontal  resolution  of  0.5° 
and  is  run  twice  daily  providing  forecast  products  at  6- 
hourly  intervals.  The  March  12-22,  1997  period  was 
chosen  to  represent  typical  early-spring  conditions 
which  included  a  storm  event  on  14  March  1997 
producing  waves  in  excess  of  4  m  as  indicated  by  buoy 
observations  and  model  results  discussed  in  Sec.  2. 1 .2, 
Figure  3  depicts  comparisons  between 
observed  winds  at  buoy  stations  41002,  41004,  and 
CMAN  station  FPSN7  (see  Figure  2  for  locations)  and 
the  NORAPS  lO-m  winds.  Figures  3a-c  show  good 
agreement  with  the  wind  direction  for  all  three 
locations  with  buoy  station  41002  showing  the  highest 
correlation  (r=0.92)  and  CMAN  station  FPSN7,  located 
at  Frying  Pan  Shoals,  with  the  lowest  value  of  r=0.82. 
The  wind  speeds  shown  in  Figures  3d-f  show 
correlations  ranging  between  0.75  and  0.85,  with  RMS 
errors  between  2.21  and  2.42  m/s.  Winds  on  12  March 
were  from  the  northeast  ahead  of  a  warm  front  which 
passed  to  the  north  of  the  area  by  9Z  on  14  March.  The 
winds  gradually  shift  from  southerly  ahead  of  an 
approaching  cold  front  and  shift  to  the  northwest  on  15 
March  near  12Z.  Winds  were  strongest  at  CMAN 
station  FPSN7  on  14  March  18Z  with  a  speed  of  20 
m/s.  Although  NORAPS  assimilates  buoy 
observations  into  its  analysis,  the  correlations  shown 
here  are  reasonable  considering  that  in  addition  to 
analyzed  fields  available  every  12  hours,  NORAPS  6- 
hr  forecasts  were  included. 

2.1.2  WAM  Results 

Typical  frequency  and  directional  energy 
distributions  from  March  1997  are  shown  in  Figures  4- 
5  from  WAM.  Figure  4  shows  the  peak  of  the  storm  on 
14  March  18Z,  with  a  peak  wave  period  of  9.6  sec,  a 
wave  height  of  4.5  m,  and  a  mean  wave  direction  of 
345®  (relative  to  north).  The  winds  at  the  peak  of  the 
storm,  that  are  generating  the  waves,  are  over  18  m/sec 
from  the  SE  and  S  (see  Figure  3).  Figure  5  shows  the 
frequency  and  directional  energy  distribution  24  hours 
after  the  storm  peak  (15  March  18Z).  The  spectral 
peak  from  the  storm  is  visible  at  a  frequency  of  0.1  Hz 
and  a  direction  of  330°,  but  a  second  spectral  peak 
around  0.17  Hz  with  a  direction  of  180®  also  appears. 
This  second  peak  corresponds  to  a  wind  event  with  a 
wind  speed  of  15  m/sec  from  the  north  (see  Figure  3). 

Figure  6  depicts  contours  of  significant  wave 
height  (m)  for  the  entire  WAM  model  domain  on  14 
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direction  and  (D-F)  wind  speed  during  March  12-22,  1997. 


Figure  4.  WAM  frequency  and  directional  distributions  of  wave 
energy  for  14  March  1997  18Z. 


March  18Z.  Maximum  heights  are  found  in  the  central 
portion  of  the  area  with  values  near  4.5  m.  Arrows 
show  the  direction  the  waves  are  moving,  generally 
from  the  SSE. 

Figures  7a-c  depict  comparisons  of  observed 
(solid  line)  versus  WAM  (open  circles)  significant 
wave  height  at  buoy  stations  41002,  41004,  and 
CMAN  station  FPSN7.  The  closest  WAM  grid  point 
to  the  buoy  was  chosen  as  the  comparison  location. 
Buoy  station  41002,  near  the  southeast  comer  of  the 
model  boundary,  had  the  lowest  correlation  at  0.90  and 
a  RMS  error  of  0.42  m.  CMAN  station  FPSN7,  closest 
to  the  STWAVE  model  grid,  showed  the  highest 
correlation  with  r  =  0.94  and  an  RMS  error  of  0.29  m. 
Although  all  three  buoys  are  outside  the  area  of  interest 
for  the  shallow-water  modeling  effort,  they- 
demonstrate  that  the  WAM  spectra  used  as  input  to 
STWAVE  and  REFDIFl  are  reasonable. 
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Figure  6.  WAM  significant  wave  height  (m)  for  14  March  1997 
18Z.  Arrows  show  direction  waves  are  propagating. 


Mean  wave  period  from  WAM  versus  observation 
is  shown  in  Figures  7d-f.  Wave  periods  range  from  4- 
10  seconds,  with  the  longest  periods  found  during  the 
peak  of  the  storm  on  14  March.  Buoy  station  41004 
and  CMAN  station  ,  FPSN7  show  reasonable 
correlations,  however  WAM  data  near  41002  are 
consistently  2  seconds  higher  than  observation.  This 
behavior  may  be  a  computational  artifact  because  buoy 
station  41002  is  near  the  southeast  comer  of  the  WAM 
boundary. 

2.2  STWAVE 

Due  to  time  step  limitations  related  to  the  grid 
resolution,  WAM  seldom  mns  at  resolutions  higher 
than  3  minutes  (i.e.  5  km).  Since  bathymetry  changes 


Figure  5.  WAM  frequency  and  directional  distributions  of  wave 
energy  for  15  March  1997  18Z. 
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more  rapidly  in  shallow  water,  a  smaller  grid  spacing 
in  necessary  for  computing  shallow  water 
The  spectral  wave  transformation  model  STWAVb 
(STeady-state  spectral  WAVE)  (Resio,  1987,  1988a, 
1988b  Davis  1992)  was  selected  to  transform  offshore 
wave  spectra  that  were  hindcast  using  WAM  mto  the 
nearshore.  This  2-D,  spectral  wave  model  was  selected 
because  it  simulates  wave  trmsformation  over  complex 
bathymetry  with  local  wind  input. 

STWAVE  numerically  solves  the  steady-state 
spectral  energy  balance  equation,  in  which  the  source 
terms  include  wind  input,  nonlmear  wave-wave 
interactions,  dissipation  within  the  wave  field,  wave- 
bottom  interactions,  and  depth-induced  breakmg.  The 
assumptions  made  in  STWAVE  are  mild  bottom 
slopes,  negligible  wave  reflection,  spatially 
homogenous  offshore  wave  conditions,  steady  wave 
and  wind  conditions,  and  linear  refraction  and  shoalmg. 

2.2.1  STWAVE  Setup 

STWAVE  operates  on  a  flat  grid  with  square  grid 
cells.  The  optimal  grid  orientation  is  for  the  y  axis  to 
be  aligned  with  the  bathymetry  contours  and  the  x  axis 
to  be  aligned  normal  to  the  contours.  This  onentation 
allows  the  greatest  range  of  offshore  wave  angles  and 

the  most  reliable  modeling  results.  rKiut  ^ 

The  Naval  Research  Laboratory  (NRL) 
supplied  bathymetry  for  the  Camp  Leje^e  region  w^ 
a  resolution  of  approximately  100  m.  The  bathymeti^ 
was  based  on  data  from  a  combination  of  sources.  (1) 
National  Ocean  Service,  (2)  high-resolution 
bathymetry  collected  in  April  1^96  by  the  Naval 
Research  Laboratory  in  support  of  the  Purple  S^ 
Exercise  (e.g.  Nichols  and  Earle  1996)  and  (3)  date 
provided  by  the  Naval  Oceanographic  Office.  A 
teiangular  irregular  network  (TIN)  was  used  to  create 
the  3  arc  sec  bathymetry.  Hurricane  Fran  ma 
landfall  just  south  of  Camp  Lejeune  in  early  September 
1996.  This  bathymetry  does  not  represent  any  local 

changes  due  to  the  hurricane.  ^ 

The  STWAVE  grid  was  generated  usmg  the 
ACES2.0  LiRGG  (Uniform  Rectilinear  Grid  GUI ) 
software  (Leenknecht  and  Tanner  1997).  Using 
URGG  the  digital  bathymetry  file  (in  the  format  of  x, 
y,  and  i,  with  x  and  y  in  geographic  coorduiates  and  z 
in  meters)  were  imported  and  the  coordmates  were 
transformed  into  Universal  Transverse  Mercator 
coordinates.  The  orientation  of  the  gnd  was  selected  to 
be  330“  relative  to  North,  so  the  grid  was  aligned  with 
the  shoreline  and  bottom  contours.  For  Camp  Lejeune 
the  grid  was  specified  as  201  cells  m  Ae  cross-shor 
and  301  cells  in  the  longshore,  with  a  grid  resolution  of 

'Graphical  User  Interface 


250  m  URGG  uses  Delauney  triangulation  to  develop 
the  grid  and  linear  interpolation  to  assip  elevations  at 
each  cell.  After  generation  of  the  grid,  the  sea  bed 
elevations  (-  values)  were  converted  to  depths  (+ 
values),  as  required  for  STWAVE  input. 

The  main  driving  for  the  nearshore  waves  is 
wave  spectra  input  on  the  offshore  bound^  of  4e 
STWAVE  grid.  These  input  spectra  are  the  output 
from  the  time-dependent  WAM  wave  model  runs,  as 
discussed  in  Section  2.1.2.  Figure  8  shows  Ae 
STWAVE  grid  and  the  location  where  the  WAM 
spectra  were  applied  at  the  offshore  boundary. 
STWAVE  is  run  with  the  same  frequency  resolution  as 
WAM  However,  the  STWAVE  grid  orientation  and 
dfrectional  resolution  differ  from  the  WAM  output.  A 
0"  wave  direction  in  WAM  is  a  wave  propagatmg  to 
the  north.  In  STWAVE,  a  0“  wave  direction  is 
propagating  normal  to  the  offshore  edge  of  die  gnd. 
Thus  the  WAM  spectra  were  translated  mto  the 
STWAVE  orientation  (STWAVE  directions  =  330“  - 
WAM  directions).  Also,  the  directional  coverage  and 
resolution  differ  between  the  WAM  output  and  the 
STWAVE  input.  WAM  spectra  cover  a  full  360  with 
a  resolution  of  15“.  STWAVE  spectra  cover  a  half 
plane  (180“)  with  a  resolution  of  5“.  WAM  spectra 
were  truncated  to  a  half  plane  (neglecting  waves 
traveling  away  from  the  coast)  and  the  resolution  was 

linearly  interpolated  from  1 5“  to  5“. 

Figure  8  shows  the  depth  contours  over  the 
STWAVE  grid  The  contours  are  in  meters,  with  0  m 
representing  the  shoreline.  The  foreshore  slope  is 
relatively  steep  out  to  a  depth  of  10  m,  where 
becomes  very  gentle.  The  contours  approxmately 
parallel  to  the  shoreline.  The  offshore  edge  of  the  gnd 
is  at  a  depth  of  30  m.  Since  the  bathymetry  at  Camp 
Lejeune  is  quite  regular  (reasonably  straight  ^ 
parallel  contours),  the  wave  transformation  is  tairiy 
uniform  along  the  shore.  The  dominant  P^o«ss“  ^ 
wave  shoaling,  refraction,  and  depth-mduced  breaking 
Wind  speed  and  direction  are  also  mput  to 
STWAVE  in  the  spectral  input  file.  The  wm 
parameters  were  supplied  from  NORAPS  (see  Section 
2.1.1).  The  NORAPS  wind  speed  and  direction  at  tne 
STWAVE  grid  offshore  boundary  (34.2®  N  and  7  ^ 
W)  were  translated  into  the  STWAVE  reference 
for  input  to  STWAVE.  The  wind  speed  and  direction 
were  assumed  constant  over  the  STWAVE  mo 

The  hydrodynamical  model  ADC1RC-2DDI 
(Luettich  et  al.  1992,  Westerink  et  al.  1994)  w^  used 
to  calculate  water  surface  elevations  which  wer  P  ^ 
into  STWAVE.  The  Advanced  ^IRCu lation  m 
(ADCIRC)  is  a  2D  finite-element  model 
tidal  constituents,  winds  and 
(from  NORAPS)  to  compute  water  surface  ele 
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ADCIRC  also  outputs  depth-averaged  currents, 
however,  they  were  not  used  in  the  simulations 
discussed  in  this  paper. 


-77.5’  -77' 


Figure  8.  STWAVE  grid  and  bathymetric  contours.  Circles 
denote  locations  of  WAM  spectra.  Inner-box  denotes  REFDIFl 
model  grid.  STWAVE  outer  boundary  is  approximately  parallel 
to  the  30m  isobath. 

2*2.2  STWAVE  Results 

The  waves  during  the  March  1997  simulation 
are  representative  of  spring  conditions  with  the  growth 
and  decay  of  waves  from  a  large  storm  and  from 
several  smaller  events.  The  wave  heights  ranged  from 
0.5  to  4.5  m  and  peak  periods  of  4  to  10  sec.  The  mean 
wave  directions  vary  throughout  the  simulation  period, 
but  in  the  largest  event  waves  are  propagating  toward 
the  shore  at  345°,  Figure  9  shows  transformed 
frequency  and  directional  distributions  of  wave  energy 
at  four  locations  on  the  STWAVE  grid  (14  March  1997 
18Z).  The  locations  represent  a  cross-shore  transect  at 
the  central  section  of  the  beach.  The  local  water  depths 
for  Figure  9  are  25.3  m,  6.5  m,  4,1  m,  and  1.6  m, 
respectively.  The  wave  directions  in  Figure  9  are 
relative  to  the  STWAVE  grid,  so  0°  represents  a  wave 
propagating  toward  330°  relative  to  north,  positive 
angles  are  more  toward  the  west  and  south,  and 
negative  angles  are  more  toward  the  north  and  east. 
The  shape  of  the  frequency  spectra  stays  about  constant 
for  each  location,  but  the  energy  level  decreases  due  to 
wave  breaking.  The  spectral  shape  stays  constant 
because  there  is  little  or  no  additional  wave  growth 
between  locations.  The  directional  distribution  of 
energy  becomes  narrower  and  more  peaked  in  shallow 
depths  due  to  refraction  turning  the  waves  so  they  are 
more  shore  normal. 

Figure  10  shows  contours  of  wave  height  for 
the  same  time  period.  The  shoreline  is  at  the  right  of 
the  plot  and  the  offshore  is  at  the  left.  The  offshore 
dashed  line  is  the  4.5-m  contour  and  the  solid  line  in 


the  center  of  the  plot  is  the  4-m  contour.  The  wave 
height  shoals  sli^tly  between  the  4.5-m  and  4-m 
contours  and  then  decays  due  to  white  capping.  Very 
nearshore,  the  wave  height  decays  due  to  depth-limited 
breaking. 


Figure  9.  STWAVE  spectra  at  25.3  m,  6.5  m,  4.1  m,  and  1.6  m 
depths  on  14  .March  1997  18Z. 


2.3  REFDIFl 

In  shallow-water,  many  dynamic  processes 
including  shoaling,  refraction,  diffraction,  and  energy 
dissipation  due  to  bottom  friction  and  depth-induced 
breaking  affect  wave  propagation.  In  the  present 
operational  Navy  suif  model,  the  RCPWAVE 
(Ebersole  1986)  wave  model  was  implemented  for 
refraction  and  diffraction  computations.  However, 
RCPWAVE  was  developed  only  for  open  coasts  with 
slowly  varying  bathymetiy.  In  some  cases,  the 
bathymetry  needs  to  be  smoothed  to  achieve  numerical 
stability.  In  addition,  it  cannot  be  applied  accurately 
for  locations  with  very  complex  bathymetry  or 
surface-piecing  features  such  as  islands,  or  semi- 
enclosed  areas.  The  coastal  wave  model  REFDIFl 
(Kirby  and  Dalrymple  1994)  has  a  more  robust 
formulation  and  does  not  suffer  from  these  limitations. 

The  REFDIFl  model  solves  the  mild  slope 
equation  with  the  parabolic  approximation  (Kirby  and 
Dalrymple,  1983).  Kirby  (1986)  extended  the  range  of 
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model  validity  for  a  wider  range  of  input  wave  angles^ 

The  model  is  solved  in  “ 

efficient  Crank-Nicholson  imphcit  scheme.  Energy 
dissipation  in  the  model  permits  treatment  of  bottom 
Snal  losses  due  to  rough,  porous  or  vtscou 
bottoms,  surface  film  and  depth-mduce  w 
breaking.  The  model  also  includes  Ae  computation  of 
wave-current  interaction,  which  is  important  at  areas 
near  inlets  and  straits. 


Figure  10.  STW.AVE  wave  height  contours  (m)  on  14  March 
1997  18Z. 

The  REFDIFl  model  is  designed  for  simulation  of 
monochromatic  and  unidirectional  wave  nam 
propagation.  For  any  realistic  wave  condit  on 
consisting  of  various 

swell  the  user  needs  to  make  mdependent  REFDIF 1 
runs  for  wave  components  at  fine  frequency  an 
direction  bandwidths.  The  results  of  these  sep^te 
runs  are  linearly  combined.  Outside  the  surf  zon 
where  depth-induced  breaking  rarely  occm,  *e 
superposition  approximation  is  .  p;  i,» 

succSsfully  applied  to  the  Southern  California  Bight 
by  O’Reilly  and  Guza  (1993). 

In  w^ve  simulation,  REFDIFl  wave  hmdcasts, 
nowcasts,  or  forecasts  can  be  used  for  each  input  wave 
condition  each  time  or  a  transfer  function  approach  cm 
be  applied.  In  the  transfer  function  app^o"*. 
calculations  are  made  for  all  possible  ^ 

angular  components,  and  the  results  are  sa 
tabular  form.  The  user  only  needs  to  derive  *e 
-function  once  for  a  given  area  assummg  no  changes  m 


local  bathymetry.  The  transfer  f^ctmn  in 

the  model  domain  consists  of  amplitude  and  phase  as 
SSrof  input  wave  conditions.  For  any  given 
wave  spectrum  input,  the  spectrum  is  first  divided  mto 
many  wave  components.  TTie  T. 

of  individual  wave  com^nents  is  modffied  by  * 
corresponding  transfer  function  and  then  Imear  y 
combLd  wi*  the  results  from  other  components  to 
derive  the  final  results. 

2.3.1  BF.FDIFl  Setup 

In  this  simulation,  the  transfer  function  approach 
i,  used.  The  model  is  run  for  every  2"  ranging  iom 
-70  to  +70  degree  relative  to  the  model  ^d. 
Frequency  covers  from  0.06  to  0.35  Hz  with  a  tota  o 
18  bands  From  0.06  to  0.2  Hz,  model  runs  are 
conducted  at  0.01  Hz  interval.  For  hi^  frequency 
waves  above  0.2  Hz,  the  frequency  mterval  mcreases  to 
0.5  Hz.  All  together,  1278  individual  wave  rans  axe 
conducted.  The  transfer  functions  are  found  to  e 
highly  variable  over  a  directional  spread  of  a  few 
degrees.  Similar  variations  have  been  reported  by 
O’Reilly  and  Guza  (1993).  Smoothing  is  applied  to  the 
transfer  functions  to  reduce  the  scatter.  In  view  of  die 
fact  that  most  of  the  energy  are  within  50“ 
and  negative)  to  the  shore  normal  for  the  time  period 

studied,  no  additional  runs  are 
arid  In  a  parabolic  model  such  as  REFDIFl,  the 
results  are  degraded  for  high  incident  angles  beyond 
50°.  In  such  cases,  grid  rotation  is  often  needed  to 
consider  high  incident  angles.  Additional  discussions 
about  REFDIFl  modeling  issues  are  covered  in  a  paper 
by  Kaihatu  et  al.  (1997)  in  this  proceeding. 

2.3,2  WF.FPIFI  Results 

The  WAM  spectra  nearest  to  the  shoreline, 
located  at  34.4°  N,  -77.2°  W,  at  a  depth  of  18.4  m,  is 
used  as  input  to  REFDIFl.  The  time  span  coven  the 
same  period  as  described  for  the  STWAVE  hmdcasts 
2  2  2)  A  sample  companson  of  frequency 
Ipectra  between  STWAVE,  ^FDIFI  (both  ^ 
depth)  and  WAM  (18-m  depth)  for  14  March  18Z  is 
shown  in  Figure  11.  The  REFDIFl  spec^  does  not 
differ  much  from  the  input  WAM  spectra.  This  is  to  be 
expected  since  waves  rarely  lose  energy  at  Ae  7-m 
depth  due  to  breaking.  The  slight  diff^t^ce  between 
STWAVE  and  REFDIFl  can  be  attributed  to  me 
inclusion  of  additional  wave  growth  and  white  cappmg 
dissipation  in  STWAVE.  The  major  chMge  on  wave 
characteristic  in  this  case  is  the  wave  angle  f /bown  m 
Figure  9.  The  evolution  of  sea  state  for  14  March  (00- 
18  Z)  is  illustrated  in  Figure  12.  The  saturate  of  hig 
frequency  waves  is  evident.  The  peak  frequency 
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changes  from  high  to  lower  frequency,  showing  a 
typical  wave  growth  due  to  wind. 


Figure  11.  Comparison  of  WAM,  STWAVE  and  REFDIFl 
wave  spectra  for  14  March  1997 18Z  near  Camp  Lejeune,  NC. 


Figure  12.  REFDIFl  spectral  density  for  March  14, 1997.  Note 
changes  in  peak  frequency  from  high  to  lower  frequency  due  to 
wind  induced  wave  growth. 

2.4  Navy  Standard  Surf  Model  /SURF96 

The  Navy  Standard  Surf  Model  (NSSM)  was 
developed  (Earle,  1988,  1989)  because  previous  Navy 
surf  forecasting  techniques  were  largely  manual,  were 
based  on  methods  dating  to  the  1950’s,  and  did  not 
adequately  consider  local  shallow-water  effects. 
NSSM  is  contained  in  the  Geophysical  Fleet  Mission 
Program  Library,  the  Tactical  Environmental  Support 
System,  and  the  Mobile  Oceanographic  Support 
System.  The  NSSM  has  been  used  extensively  by  the 
Fleet  and  for  several  applications  such  as  estimating 
climatological  surf  conditions  at  selected  beaches  for 
development  of  the  Marine  Corps’  Advanced 
Amphibious  Assault  Vehicle  (McDermid  et  al.  1997; 
Nichols  et  al.  1997) 

The  NSSM  is  composed  of  two  main  modules,  a 
wave  refraction  model  and  the  surf  model  itself, 


SURF96.  The  most  recent  version  of  the  surf  model, 
SURF96  incorporates  several  theoretical  and  numerical 
improvements  (Earle  et  al.l997;  Hsu  et  al.  1997)  and 
was  used  for  this  work.  The  surf  model’s  strengths  are; 
rapid  operation  for  field  use  and  simulations,  relatively 
simple  operation  and  mathematical  robustness  for  non¬ 
expert  use,  realistic  depiction  of  breaking  wave 
locations,  and  provision  of  detailed  information 
(breaker  heights,  types,  longshore  currents)  across  the 
surf  zone.  The  model  automatically  provides 
parameters  for  planing  amphibious  operations  as 
described  in  the  Joint  Surf  Manual 
(COMNAVSURFPAC/COMNAVSURFLANT,  1987). 

SURF96  is  a  one-dimensional  parametric  model 
largely  based  on  concepts  developed  by  Thornton  and 
Guza  (1983,  1986).  The  assumptions  include:  (1) 
approximately  straight  and  parallel  bottom  contours 
within  the  surf  zone;  (2)  a  directional  wave  spectrum 
that  is  narrow  banded  in  frequency  and  direction;  (3)  a 
Rayleigh  wave  height  distribution;  and  (4)  linear  wave 
theory.  The  directional  wave  spectrum  used  for  model 
initialization  can  be  obtained  from  a  wave  model 
(WAM,  STWAVE,  REFDIF)  or  measurements.  If  an 
offshore  directional  wave  spectrum  is  used,  it  must  be 
properly  refracted,  shoaled  and  diffracted  to  a  starting 
depth  which  lies  seaward  of  the  breaker  line.  The 
directional  wave  energy  distribution  js  reduced  to  three 
representative  physical  values;  (l)~ffle  direction  of  the 
vertically  averaged  wave  momentum  flux,  (2)  the 
incident  wave  energy,  and  (3)  the  dominant  wave 
frequency.  The  direction  of  the  wave  momentum  flux 
is  computed  using  radiation  stress  (e.g.  Longuet- 
Higgins,  1970a,b).  The  model  incrementally  calculates 
wave  energy  from  which  wave  heights  are  determined 
along  a  transect  normal  to  the  beach  to  very  near  the 
beach  and  still  water  level.  As  waves  move  through 
the  surf  zone,  the  average  rate  of  energy  dissipation 
due  to  wave  breaking  and  frictional  dissipation 
balances  the  gradient  of  the  shoreward  energy  flux. 
Longshore  current  calculations  at  each  increment  are 
based  on  longshore  current  theory  using  radiation  stress 
(e.g.  Longuet-Higgins,  1970a,b).  Current  calculations 
include  local  wind  stress  effects.  The  longshore 
current  module  is  being  examined  to  provide  improved 
performance  for  beaches  with  shallow  offshore  bars. 
Such  a  bar  was  not  present  for  these  data  (Figure  13). 

The  percentage  of  breaking  and  broken  waves 
across  the  surf  zone  is  a  model  output.  Surf  zone  width 
is  considered  as  the  most  seaward  point  where  ten 
percent  of  the  waves  are  breaking  or  broken. 
Percentages  of  breaker  type  (spilling,  plunging,  or 
surging)  are  determined  across  the  surf  zone  from  the 
probability  distribution  of  breaking  waves  using  a 
widely  accepted  parameterization  based  on  breaker 
height,  breaker  period,  and  bottom  slope.  Breaker 
height  parameters  specified  in  the  Joint  Surf  Manual 
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are  determined  for  regions  where  breakers  are  largest. 
The  Modified  Surf  Index  (MSI)  defined  in  the  Joint 
Surf  Manual  is  calculated  automatically.  The  MSI 
determines  whether  particular  types  of  conventional 
landing  craft  (not  air-cushion  vehicles)  can  be  used. 

2.4.1  SURF96  Results 

SIJRF96  hindcasts  for  the  simulation  period 
12-22  March  1997  were  performed  using  wave  spectra 
from  two  sources:  (1)  REFDIFl  and  (2)  STWAVE.  A 
Navy  SEAL  beach  survey  collected  for  the  Purple  Star 
Exercise  was  used  as  the  beach  profile  in  both  sets  of 
SURF96  model  runs.  Hourly  wind  speed  and  direction 
fromlhe  Marine  Corps  Air  Station  (34.7°  N,  77.43°  W) 
at  New  River,  NC  were  used.  NORAPS  winds  were 
available,  however  it  was  felt  that  observed  winds 
would  provide  the  most  realistic  input  into  SURF96. 
Because  water  elevations,  including  tides,  affect  surf, 
water  elevations  from  ADCIRC  were  input  into 
SURF96.  The  SURF96  code  was  adapted  to  accept 
wave  spectra  from  REFDIFl  and  STWAVE. 

Shallow-water  wave  spectra  from  STWAVE 
and  REFDIFl  hindcasts  are  saved  for  selected 
nearshore  points.  In  this  paper,  the  directional  spectra 
from  both  models  at  a  7-m  depth  along  the  beach 
survey  line  are  used  as  input  to  SURF96.  The  7-m 
location  is  selected  so  that  only  a  small  probability  of 
depth-induced  wave  breaking  will  occur  at  this 
location.  This  location  was  selected  to  assure  that  the 
assumption  of  no  depth-induced  breaking  is  valid  for 
the  REFDIFl  transfer  function  approach. 

Figure  13  depicts  SU]^96  results  for  15 
March  03Z.  The  top  panel  shows  the  Navy  SEAL 
beach  team  survey  used  as  input  to  SURF96.  The  slope 
is  rather  steep  from  the  coastline  to  approximately  250 
meters  offshore  (at  a  depth  of  4  m)  where  the  slope 
levels  off.  Figure  13b  shows  significant  wave  height 
along  the  shore-normal  transect  to  900  m  offshore.  The 
solid  line  depicts  SURF96  results  with  input  from 
STWAVE  and  the  open  circles  represent  results  in 
which  REFDIFl  spectra  were  input.  Figure  13c  shows 
almost  identical  results  of  the  percentage  of  breaking 
waves  when  comparing  SURF96  output  using  inputs 
from  REFDIFl  and  STWAVE.  The  surf  zone  width 
(not  shown)  is  defined  as  the  offshore  location  where  at 
least  10%  of  the  waves  break.  This  occurs  at  a  distance 
of  approximately  225  m.  (Figure  14  depicts  the 
variation  of  surf  zone  width  with  time).  The  waves  do 
not  begin  to  break  until  water  depths  become  less  than 
4  m  because  the  input  significant  wave  heights  are 
between  2  m  and  3  m.  Finally,  Figure  13d  shows  the 
longshore  current  variation  with  distance  and  model 
input  into  SURF96.  The  maximum  longshore  current 
occurs  at  a  distance  of  150  m  offshore.  Currents  using 
STWAVE  input  show  larger  magnitudes  than  those 


using  REFDIFl  input  in  Figure  13  because  STWAVE 
surf  zone  incident  wave  angles  are  slightly  larger  than 
those  for  REFDIFl.  Wave  directions  are  incident  at 
small  angles  from  a  peipendicular  to  the  beach  so  that 
currents  are  sensitive  to  small  changes  in  wave 
directions  just  outside  of  the  surf  zone.  At  various 
times,  slightly  different  wave  directions  rather  than 
different  wave  heights  seem  to  account  for  current 
differences  between  the  two  model  inputs. 

Figure  14  presents  time  series  of  SURF96 
output  during  12-22  March  1997.  The  top  panel  depicts 
wind  speed  and  direction  provided  from  the  Marine 
Corps  Air  Station.  Wind  speeds  are  highest  on  14 
March  in  advance  of  a  strong  cold  front  as  wind 
directions  are  generally  from  the  south.  The  wind 
gradually  veers  to  the  west  and  northwest  as  the  cold 
front  passes  through  the  area  and  become  southerly  on 
18  March.  Late  in  the  day  on  19  March  the  winds 
again  shift  to  a  direction  from  the  north  as  yet  another 
cold  front  passes  through  the  area. 

TTie  second  panel  depicts  the  significant  wave 
height,  which  is  an  input  parameter  to  the  surf  model. 
The  solid  line  represents  the  STWAVE  whereas  the 
open  circles  represent  REFDIFl  output.  Both  models 
yield  very  similar  results.  The  highest  wave  heights 
are  associated  with  the  major  storm  event  occurring  on 
14  March. 

The  third  panels  shows  the  wave  angle  at  the 
7-m  depth.  0°  represents  waves  moving  perpendicular 
to  the  shoreline.  A  positive  breaker  angle  is  shown  to 
be  for  breakers  moving  toward  the  right  flank  relative 
to  a  sight  line  perpendicular  to  the  beach.  Again, 
results  from  SURF96  with  REFDIFl  and  STWAVE  as 
inputs  are  similar.  From  12  March  21Z  to  13  March 
18Z,  REFDIFl  shows  a  wave  direction  approximately 
10°  more  negative  than  STWAVE.  Some  of  this 
variation  may  be  attributed  to  wind  forcing,  which  is 
included  in  STWAVE,  but  is  not  for  REFDIFl. 

The  fourth  panel  depicts  the  surf  zone  width, 
which  is  defined  as  the  offshore  location  where  at  least 
10%  of  the  waves  are  breaking.  Throughout  the  10- 
day  hindcast  period,  surf  zone  widths  are  generally 
between  60-200  m.  The  exception  occurs  during  an 
18-hr  period  from  14  March  15Z  to  15  March  03Z 
when  the  width  exceeds  300  meters,  when  a  major 
storm  event  affected  the  area.  The  effects  of  tides  are 
evident  by  the  modulation  in  surf  zone  width 
throughout  the  period. 

The  fifth  panel  shows  the  longshore  current 
with  magnitude  maxima  near  0.8  m/s  on  14  March. 
There  is  a  change  in  current  direction  associated  with 
changes  in  wind  and  wave  directions  relative  to  the 
beach. 

The  final  panel  shows  the  modified  surf  index 
(MSI)  which  is  derived  using  the  maximum  wave 
height,  maximum  longshore  current,  breaker  angle, 
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*  A  onH  breaker  tvpe  and  based  on  the  Joint 

::r,:;^“'Msror of 

conventional  amphibious  landmg  craft. 

3,0  SZER  Model  Validation 

All  the  models  contained  in  the  SZER  model 
suite  have  been  validated  and  tested  in  either  laboratory 

:;rrfreldtests^  a;^197T^ ^^r/a^ 

SllSVouck,  NC  during  the  Duck  Experiment  on 
Low-Frequency  'and  uSd 

S  rco“;trs”uen  ..s.. . ; 

“"^%SX'p”riotof26Aug„s.-5Sep»j^^^^ 
1997  surf  pone  model  validation  data  were  collected 
Sw  during  the  Joint  Countermeasures 

hs  Bri:r«="w- 

'buoys  developed  by  the  NRL  j 

OrpanoffraDhv  Warfare  Support  (TO  ) 

Ote  Space  and  Naval  W^ate  ^To  m  oTwamt. 
(SPAWAR)  were  moored  m  3,  7  lU  m 
In  addition,  pressure  sensor  wave  and  ti  g 

Hpnloved  outside  and  withm  the  surf  zone,  me  sun 
deployeo  ouisi  .jpo  ^ere  monitored  usmg 

zone  environmental  processes  w  «iinw 

X  image  processing  techniques 

quantification  of  both  ® ^ 

Jaratnetets  including  wave  period,  wave  »g  a. 

Leed  surf  zone  width,  shoreline  location,  and  sand 
rt™.  Estimates  of  '""Ssho^  cPtTm  an 
bathymetry  can  also  be 

plans  call  for  running  „ea  and 

described  in  this  paper  for  ..^^Lh  during  this 

validating  model  output  with  data  collected  durmg 

exercise. 

4.0  SUMMARY 

A  coupled  wave/surf  model  suite  consisting  of 
“off-the-shelf’  v/ave  and  surf  models  is  evaluate  . 
series  of  hindcasts  are  performed  for  e  n 
NC  area  during  the  period  12-22  March  1997. 

The  deepwater  WAM  model  shows  ry 
good  agreement  wL  available  buoy  and  CMAN  dato^ 
WAM  directional  wave  spectra  are  J?P 

into  two  shallow-water  wave  models:  (1)  STWA 
Pand  (2)  REFDIFl.  Model  outputs  from  ST  WAVE  an 
;  REFDIFI  are  remarkably  similar. 


Outputs  from  STWAVE  snd  REFDlFi  dttve  a 

zecent  vSsten,  SURFW,  of  the  NSSM  to  prov,^ 

detaiied  information  across  *e  “nf  military 

♦oi-c  nf  interest  for  simulations  of  military 

‘’T  Id  ajSfr ious  operations  are  similar  for 

systems  and  mphib'ous  p^^^  elevations, 

inputs  from  AnriRC  illustrates  effects  that 

width.  Relativ  y  ^ 

widths,  and  staulate  conditions 

Zc-rm"-  of  miiitaty  systems  and 
amphibious  operations. 
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